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ABSTRACT 

Results are presented for a preliminary experiment involving satur- 
ated pool boiling at 1 atmosphere from rotating spheres, 2 and 3 in. in 
diameter, which were immersed in liquid nitrogen (LN 2 ). Additional 
results are presented for a stationary, 2-in. -diameter sphere, quenched 
in LN 2 , which have been obtained utilizing a more versatile and com- 
plete experimental apparatus that will eventually be used for addi- 
tional rotating sphere experiments. The speed for the rotational tests 
was varied from 0 to 10,000 rpm. The stationary experiments parametri- 
cally varied pressure and subcooling levels from 0 to 600 psig and from 
0 to 50°F, respectively. During the rotational tests, a high speed 
photographic analysis was undertaken to measure the thickness of the 
vapor film surrounding the sphere. The average Nusselt number over the 
cooling period was plotted against the rotational Reynolds number, sta- 
tionary sphere results included local boiling heat transfer coefficients 

at different latitudinal locations, for various pressure and subcooling 
levels. 


INTRODUCTION 

The main shaft ball bearings in the Space Shuttle Main Engine (SSME) 
high pressure oxidizer turbopumps (HPOTP) are exposed to high tempera- 
tures because of frictional heating during operation, resulting in bear- 
ing ball discoloration. To improve the heat removal mechanism from these 
balls, it is necessary to better understand the boiling heat transfer 
phenomenon associated with rotating balls in a cryogenic fluid. Further- 
more, an exhaustive literature search was conducted to identify heat 
transfer correlations related to pool and forced-convection boiling from 
quenched spheres. Several correlations were found, but none encompass 
all of the independent parameters required for this application such as 
rotational speed and pressure. As a result, an extensive experimental 
investigation was initiated to better understand the effect of rotation, 
coolant pressure, and subcooling level on boiling heat transfer from 
spherical surfaces. 

The primary objective of this experimental program was to parametri- 
cally study the effect of ball rotational speed, coolant pressure, and 
subcooling on boiling heat transfer from a heated sphere quenched in a 


*Work reported herein was sponsored by NASA/Marshall Space Flight Center 
under Contract NAS8-40000. 
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liquid nitrogen (LN2) bath. The experimental program used spheres of 
two different sizes to obtain ball surface velocities and centrifugal 
forces in the same range as those experienced in the HPOTP ball bearing. 
In addition, the ball was made of 440C CRES material, which is identical 
to the engine's bearing balls. Because of safety reasons, liquid nitro- 
gen was used as the test liquid, as the thermodynamic and transport 
properties are similar to those of liquid oxygen (engine coolant). The 
range of pressures studied was varied from atmospheric to supercritical. 
The effect of cross flow was not included in this first phase of the 
study, although it is recognized that it eventually must be incorporated 
into the experiments in order to closely simulate the SSME bearing 
environment. The rotational speed was varied from 0 to 10,000 rpm, at 
which speed the 3— in. ball rotated at 67 percent of the equivalent 
surface speed for the 0.5-in. engine ball. 

TEST APPARATUS 

The early quench tests were carried out with a rotating sphere 
under atmospheric conditions. Since thermocouples were not attached, 
only the total elapsed quench times until the ball cooled to the point 
of transition boiling were recorded. The second phase of testing under 
pressure required a much more involved test facility. This facility is 
capable of testing over a range of pressures, rotational speeds and 
levels of liquid subcooling. Both of these facilities were designed 
and built for quench tests in a pool of LN 2 - 

First Phase 


The first phase study basically provided overall cooling times for 
a rotating sphere quenched in LN 2 - The test apparatus (Fig. 1) for 
this phase consisted of a dewar, motor, and ball. The glass dewar was 
a 10— in. diameter , 9— gal, double-wall, vacuum-sealed flask. Quenching was 
accomplished after heating the sphere in air and then raising the LN 2 
glass dewar to immerse the sphere. A hydraulic unit was used for this 
purpose. 440C stainless steel spheres (2- and 3-in. diameter) were used 
in these tests. A 3/8— in. OD., hollow. Inconel 625 shaft connected the 
motor to the ball. A 1/2 hp variable speed motor, capable of achieving 
speeds up to 10,000 rpm, was used to rotate the shaft and the sphere. 

High speed motion picture cameras were also employed in this phase 
of testing to observe the physical phenomenon along with the quantita- 
tive measurements. The photographic analysis was accomplished using 
both Fastax and Hycam high speed motion picture cameras to record the 
boiling phenomenon. Frame rates from 400 to 10,000 frames/sec were 
used in these experiments. Various combinations of front and silhouette 
back lighting were tried along with various fields of view to optimize 
the photographic results. 

Second Phase 


The second phase test facility was considerably more complex than 
that of the first phase. The liquid container had to maintain super- 
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Fig. 1. Rotating Sphere Test Facility 


critical pressures in addition to allowing for rapid pressurization 
after the Sphere was heated. Furthermore, the pressure vessel had to 
be raised and lowered as in the first phase, while the pressure in 
container was maintained relatively constant. Also, a pressurize , 
low- flow, liquid-feed system connected to the inner container was 
needed to maintain a constant pool temperature during sphere cooldown 
Finally, windows in the container were required for both test monitor g 
a^ photographic analysis. Figure 2 shows an overall view of the 
apparatus employed in this study. 

A pressurized, cylindrical steel container (10-in. diameter) was 
specifically made for this test program. An open-glass beaker, ' 7 

in diameter and having a capacity of 6.5 gal, was placed _ ^ 

oressure vessel to hold the liquid. The pressurized feed system was 
used both for filling and for maintaining a constant bath temperature. 
The slow-moving coolant was continuously introduced into the bottom o 
JE gSs container, and was allowed to overflow into an annulus formed 
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Fig. 2. Overall View of Stationary Sphere Test Facility 


between it and the pressure vessel inside wall. The surplus LN 2 was 
then evacuated from the vessel base to a receiving tank. The feed 
system was capable of operating at all pressures and liquid tempera- 
tures. Liquid nitrogen was furnished to the test supply tank, with a 
capacity of 33 gal, from an external storage tank. The test supply 
tank and feed lines were insulated. The LN 2 temperature in the 
supply tank was controlled both by pressure and the amount of warm 
gaseous nitrogen (GN 2 ) injected into the bottom of the tank. The 
temperature of the LN 2 slowly rose after pressurization, advancing 
toward the equivalent saturation temperature. Figure 3 shows the 
schematic of the flow system. A highly polished, 2-in. diameter sphere 
made of 440C stainless steel material was used for the boiling experi- 
ments. Sheathed and grounded thermocouples were used to measure the 
ball temperature-time history during cooldown as well as to estimate 
the heat transfer from the surface of the ball. The thermocouples were 
radially positioned in pairs at various latitudinal locations on the 
sphere. One thermocouple from each pair was located flush with the 
surface and the other was located radially inward between 0.20 to 0.5 
in. depending on the location. The sphere was X-rayed to accurately 
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Fig. 3. Schematic of Stationary Sphere Test Rig 
With LN 2 Feed System 

determine the thermocouple positions. Nine thermocouples, 0.02-in. in 
diameter 6 were brazed into holes made in the sphere. The thermocouple 
wires were routed through the hollow. 3/8-in. diameter. Inconel 
625 shaft to the data acquisition system. When the rotational test 
phase begins, the thermocouple will be attached to slip ring connectors, 
which in turn will be connected to the data acquisition system. Figure 
4 shows the sphere instrumentation locations. 

In addition, three thermocouples were placed in the liquid bath at 
elevations above, below, and at the sphere’s equator to check the ln 2 
bulk temperature in the vicinity of the sphere. Two thermocouples were 
also inTalled inside the supply tank to control the subcooling evel 
of the nitrogen before each test. Pressure transducers were placed 
inside the pressure vessel and in the supply tank. 

The data acquisition system consisted of two major components: (1) 

the reference Junction box and (2) the Hewlett Packard system 4000 
acauisition system. The thermocouples leads from the sphere wer 
connected to the reference junction box to provide a common junction 
for accurate temperature readings. The reference temperature was 150 
+ 1° F . The output of the box, in millivolts, was in turn sent to the 
system 4000. The scanning rate of the system 4000 was 30 channels/sec. 
The system 4000 recorded the 18 test parameters (temperature and pres- 
sure measurements) at every 1.97-sec interval. Tnternallyktheinput 
signal passed through three stages. First, the signa ^ "as ^ 

a factor of 10 using an amplifier. The precalibration amplifier s accur 
acy is ±6.1 microvolts, which corresponds to approximately ±0.3 F 
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Fig. 4. Sphere Instrumentation Locations for 
Stationary Testing 

during most of test duration, and to +0.6°F in the nonlinear domain 
which occurred briefly near the end of the test. After the amplifica- 
tion, the signal was digitized using a 14 bit analog to digital con- 
verter. Finally, the digitized data was transferred into a personal 
computer's floppy disk and eventually transformed into an IBM format 
ASCII file for posttest processing. 

To minimize the overall uncertainly of the data acquisition system, 
a customized calibration was performed with a temperature range varying 
from -320 to 212°F. A Doric unit with type-T thermocouples was the 
standard utilized for the calibration. In addition to four fixed 
temperature settings: saturated liquid nitrogen (-320°F), a mixture 

of dry ice and alcohol (-100°F), ice water (-32°F), and boiling water 
(212°F) , intermediate steps between -320°F and -232°F were traversed. 
The technique consisted of pressurizing the liquid nitrogen from satur- 
ation to its critical pressure in increments of 50 psia, letting it 
reach the equivalent saturation temperature at each step. Calibration 
data revealed that thermocouples temperature readings agree closely 
with the above-mentioned standard used and hence only minor adjustments 
were required. Following these adjustments, the accuracy of the thermo- 
couple improved to approximately +0.15°F. 
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TEST PROCEDURE 


Rotating Tests 

The glass dewar was filled with L» 2 such that the liquid level was 
at le.st\ in. above the sphere The Initial -P-et^ of. ^ sphere 

This termination time was rather easy to recogn , increased 

uous vapor flu was In "»g 

M^Prlo^to' thls S tlme (during film boiling), the bulh 

b^ ITLnd" polished^°as deeded and carefully cleaned with either 
freon, alcohol, or acetone. 

FILM BOILING REGION 



90 5 SECOND 


90 SECOND 


TRANSITION FROM FILM TO NUCLEATE BOILING REGION 



Fig. 5. Flow Visualization Results From 2-in Rotating Sphere 
(Rotational Speed - 6500 rpm) 


Stationary and Pressure T ests 

The pretest preparation was rather complex for 1 J is Pressure 

tem was pre-seroed and the hydraulic lifting unit turned ™! C P "L 

SZ 
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flow into the bottom of the inner glass container. Recording of the 
data was initiated prior to start of the flow and maintained until the 
end of the test, at which time the sphere surface temperatures reached 
the liquid temperature. The liquid nitrogen subcooling level and pres- 
sure were independently and steadily controlled throughout the test. A 
VHS video system allowed remote viewing of the sphere through one of the 
three glass viewing windows. A videotape was produced to record each 
experiment in its entire duration. 

RESULTS AND DISCUSSION 


Rotating Tests 

Over 80 trials involving 2- and 3-in. spheres were conducted, vary- 
ing the rotational speed from the stationary condition up to 10,000 rpm, 
in 1000 rpm increments. Thirty-seven of those trials were used for data 
reported herein. The remainder were used for repeatability, verifica- 
tion, photographic coverage, and apparatus checkout. 

A photographic analysis was undertaken with both Fastax and Hycam 
high speed motion picture cameras with a variety of rotational speeds, 
camera frame rates, and different fields of view. An analysis of the 
film was used to determine the average thickness of the wavy vapor film, 
encompassing the 2-in. sphere. This analysis showed that the thickness 
is between 0.030 and 0.050 of an inch, for the 6000 rpm rotational 
speed. Furthermore, it was observed that the stable film boiling 
collapse always occurred at the top of the sphere first, and quickly 
progressed downward over the sphere. 

The average heat transfer rates were computed by considering the 
sphere as a lumped-parameter system since the sphere experienced a 
relatively mild thermal transient during film boiling, which constituted 
the bulk of the test time. 

The lumped-mass equation was used to find the average film heat 
transfer coefficient (h) as a function of the beginning and final 
sphere temperatures, and the cooling time. The end temperature was 
assumed to correspond to the minimun heat flux, as once the continuous 
vapor film is disrupted, the fluid friction between the rotating sphere 
and the LN 2 increases significantly causing the liquid in the dewar 
to rotate rapidly. The sphere's temperature at this minimum heat flux 
can only be assumed to be the same as that for a stationary sphere, 
which was recorded to be approximately 70°F above the saturation tem- 
perature. Furthermore, Weaver (Ref. 1), did show that there was no 
significant shift in this temperature, corresponding to the minimum 
heat flux, with rotational speed for a sphere inmersed in freon 113. 

Calculated values of h were made dimensionless and the resulting 
average Nusselt number values were plotted against the rotational 
Reynolds number, which was based on the sphere's tangential velocity 
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and nitrogen’s vagor kinematic viscosity. The following equations were 
used to calculate h , Nu and Re 


pVC T - T 

E i n -1 a- 

AAt T o -T s 


(where p and C p are based on 440C CRES material properties.) 


v k 


R U D 


The data are shown in Fig. 6 . As noted, the Nusselt number appears 
to be insensitive to the Reynolds number up to a value of approx ima e y 
°6 men the Nusselt number begins to rise as the Reynolds number 
increases This is felt to occur when the rotational velocity or iner 
effect begins to be important with respect to the buoyancy force as 
heat removal mechanism. 



Fig. 6. Experimental Results of Rotating Noninstrumented 
2- and 3-in. Spheres 
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The determination of the value of this critical Reynolds number as 
well as the trend of the data in the rotationally dominated regime is 
considered below. 

Kreith (Ref. 2) has confirmed that for (Gr/Re 2 ) < 1/10 , the natural 
convection effects can be ignored for single-phase heat transfer from a 
rotating sphere. This same criterion is used here for the pool film 
boiling (two-phase flow), where both the inertia and bouyancy effects are 
based on vapor properties. The use of vapor properties is based on 
observations of the liquid nitrogen bath, appearing nearly motionless, 
when a rotating sphere surrounded by a vapor film, is rotating within 
it. The following expression defines the dimensionless ratio 
(Gr/Re 2 ) v : 


g <Pf - P v ) D 3 



which can be rewritten as: 


P l 



gD 


( 1 ) 


( 2 ) 


where the following approximation was made above: ( P $ - P v ) s Pf - 

Figure 6 shows the values of Re v based upon Eq. (2) and the dimen- 
sionless ratio (Gr/Re2) v < 1/10 . When this inequality is satisfied, the 
bouyancy effects are unimportant for both the 2— and 3— in. diameter 
spheres. As can be seen, the predicted values of Re v , where this 

transition occurs, agree quite well with the data. 

It should be noted that earlier work with flow film boiling over a 
stationary cylinder, as discussed by Hasan et al. (Ref. 3), the criter- 
ion Uoo/tgD) 1 ^ 2 > 10 was used instead. This criterion results from the 
assumption that the dimensionless ratio (Gr/Re 2 ) is based on the liquid 
properties for P and V . As seen from Eq. (2), this would imply that the 
equivalent (Gr/Re 2 ), based on vapor properties, be equal to (Pf/P v )/100, 
which for the case of LN 2 at atmospheric condition leads to (Gr/Re 2 ) v - 
1.5. This value is much greater than 0.1 used in Fig. 6, indicating 
that the critical Reynolds number, at which point the rotationally 
dominated regime begins, is underestimated by an order of magnitude 
using Hasan et al.'s criterion. 
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Comparison With Existing Correlations 


The only possible comparison of the present data with existing 
correlations can be made by replacing the free stream velocity (Uoo) 
with the tangential velocity (U t ) in a forced flow boiling correla- 
tion. Two separate equations for saturated flow film boiling over a 
stationary sphere are included below. 

Witte (Ref. 4) developed one of those equations from theory and 
modified the constant with test data. The resulting equation is 



(3) 


Witte’s theoretical equation had the constant 0.698 instead of 2.98. 

Another equation obtained from analysis by Epstein and Hauser (Ref. 
5) is based upon the stagnation region, but is used here for comparison. 
Their equation is 



Both equations are plotted on Fig. 6, where the Reynolds number is based 
on U t , rather than Uqq . Surprisingly, the two equations bracket the 
rotating sphere data, which lends some credibility to the substitution 
of U t for Uqo in the the Reynolds number, as a method of finding the 
film coefficient for a rotating sphere using a conventional forced con- 
vection correlation. However, it should be noted that, there is no 
solid basis for replacing the free stream velocity in a flow boiling 
correlation with the tangential velocity for a rotating sphere in a 
stagnant pool. Furthermore, it should be reiterated that, for the film 
coefficient values reported here, it was assumed that h remained con- 
stant over the entire cooling period in the film boiling regime. These 
assumptions may contribute to deviations observed between the correla- 
tions and the test data, although fairly good agreement was noticed. It 
should be noted that Irving and Westwater (Ref. 6) showed that boiling 
results were quench-rate independent for spheres 2 in. in diameter or 
larger . 
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STATIONARY AND PRESSURE TESTS 


In excess of 70 trials were conducted using an instrumented, 2-in 
sphere and varying the liquid nitrogen pressure, subcooling level, and 
sphere’s initial temperature. The pressure was varied from 0 to 600 
psig (supercritical), the subcooling from 0°F (saturation) to 50°F, and 
the sphere’s initial temperature from 70°F (room temperature) to 450°F. 
Table 1 shows the test matrix with selected tests accomplished for this 
phase . 


Table I. Text Matrix for Stationary Sphere 


Subcool i ng 







Temperature 


Pressure (psig) 




(F) 

0 50 J Too 

[ 200 

j 300 | 

400 

500* 

600* 

0 

1-5 1-62 1-18 

1-1 1 

1-46 

1-44 




1-63 1-47 

1-49 

1-56 

1-50 




1-67 



1-75 



10 

1-61 1-25 

1-19 

1-26 

1-1 / 

1-64 

1-35 


1-43 

1-30 

1-41 

1-33 


1-60 


1-53 

1-34 

1-58 

1-51 



20 


1-4 

1-24 

1-12 


1-15 



1-29 

1-31 

1-40 





1-42 


1-59 



30 


1-32 

1-23 

1-21 

1-65 

1-14 




1-36 

1-27 


1-16 





1-5/ 



40 



1-10 

1-13 


1-28 




1-22 

1-3/ 






1-55 

1-39 







1-45 







1-48 







1-52 



50 




1-6 

1-7 

1-66 





1-38 



*Subcoo 1 i ng 

level based on the saturation 

temperature at 



| critical pressure 







Local film coefficients were determined based upon the radial tem- 
perature gradient, measured from several pairs of thermocouples located 
at various latitudinal locations. The data reduction was based on an 
energy balance applied to an element encompassing each thermocouple 
pair. The expression for the heat flux per unit area is given as: 


q/A 


k 

i 


TC , I 


- T. 


TC 




C H AT, 


TC,S 

At 


(5) 


where the thermodynamic properties for the sphere are based on 440C 
CRES material for k, C p and p. 
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Figure 7 displays typical test data in terms of film coefficient 
(h) versus sphere* surface temperature for the 

eauator, for various pressure levels. It can be seen that film coeffi 
cient increases as the LN 2 pressure is increased from saturation o 
supercritical levels. The major changes in h occur when the pressure 
increases from 0 to 200 psig, and subsequently the rate of c ^nge 
less as the critical pressure is reached. Also, there is a sharp 
crease in h as the pressure becomes supercritical. 
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Fiq. 7. Film Coefficient vs Surface Temperature for 
Various Pressure Levels at 10°F Subcooling 
For a Location 30-deg Below the Equator 

Figure 8 shows test date expressed In the form of ’ “ , £ ,“? C cSfl- 
wall temperature for various LU 2 subcooling levels. The film =°« ££1 
dent Is modestly affected by changes in subcooling level; although 
there appears to be a negligible change above the 30 degrees subcooled 

condition. 

Comparison of Stationary Sphe re D ata With Poo l 
Film Boiling Correlations 

Few correlations associated with subcooled pool film boiling have 
been found in the literature. Two of those found, based upon boiling 
from a stationary cylinder, are compared with the present data for a 
sphere. Witte (Ref. 4) discussed similarities In film boiling for both 

geometries. 
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Wall Temperature (°F) 

Fig. 8. Film Coefficient vs Surface Temperature at Constant 
Pressure for Various Subcooling Levels at 30 Degrees 
Below the Equator (300 psig) 

The correlation given by Sakurai, et al (Ref. 7) was examined and 
found to show a decrease in the film coefficient with an increase in 
the subcooling level. This contradicts other data found in the liter- 
ature as well as data presented herein. Consequently, no other con- 
sideration was given to this correlation. 

On the other hand, Siviour and Ede's correlation (Ref. 8) was found 
to predict trends shown by other data as well as that of this study. 
Their equation, neglecting radiation losses is: 

v °- 613 sr AT » ( Gr v "v) 1 " 4 ( fe) 

( 6 ) 

k f / 2\l/4 

* °- 59 r AT sc ( or « et i / 

where the terms in the above equation are defined in the Nomenclature 
section. 

Results obtained with Eq. (6) are shown as dashed lines in Fig. 7 
for pressure levels of 100, 200, and 300 psig and with a subcooling 
level of 10°F. However, both constants in Eq. (6) were adjusted 
upward to fit the test data. The correction factor used is equal to 
1.6. Considering the difference in geometry between a cylinder and a 
sphere, such an adjustment would not be surprising and, in fact, Witte 
(Ref. 4) utilized it. As noted in Fig. 7, the effect of pressure on h 


155 



from results obtained with the modified Eq. (6) is in reasonably good 
agreement with the trend displayed by the test data. However, the test 
data begin to deviate significantly from the correlation when the wall 
temperature is below -100°F. 


CONCLUSIONS 


The rotating sphere saturated pool film boiling test results tend 
to agree with results obtained from forced flow boiling correlations 
when tangential velocity replaces crossflow velocity. 

The subcooled pool film boiling test data for a sphere show a sub- 
stantial increase in the film coefficient with pressure rise through 
supercritical pressures . 

The subcooled pool film boiling test data for a sphere exhibit a 
modest increase in the film coefficient as the subcooling level 
increases . 


NOMENCLATURE 


A 

c p 

D 

g 

Or 

h 

*fg 

hfg 

k 

k 

ln 2 

NU 

Pr 

R 

Re 

Sp 

Tb 

Tf 

To 

T S 

T TC,1 

t tc,s 

T W 

At 

ATop 

at tc , s 

AT W 

U 

V 


area of sphere 

specific heat 

diameter of sphere 

gravitational constant 

Grashof number 

film coefficient 

latent heat of vaporization 

hfg — hfg + 0.5 Cp f y/ ATyj 

thermal conductivity 

thermocoup le spac ing 

liquid nitrogen 

Nusselt number 

Prandtl number 

radius 

Reynolds number 

sp - c p ,v AT SC /h fg 
bulk temperature 
final sphere temperature 
initial sphere temperature 
saturation temperature 

thermocouple temperature inside the sphere 
thermocouple temperature at the surface 
wall temperature 
elapsed time 
T s - T b 

surface thermocouple temperature drop over time period At 

T w - T s 

velocity 

volume of sphere 
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Subscripts 


i - liquid 

t =» tangential 

v - vapor 

Greek Symbols 

» kinetic viscosity 
- rotational speed 
« density 
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